hep-th/0609100 



Brane Cosmological Evolution with Bulk Matter 



C. Bogdanos and K. Tamvakis 

Physics Department, University of Ioannina 
Ioannina GR451 10, Greece 

Abstract 

We consider the cosmological evolution of a brane for general bulk matter content. In 
our setup the bulk pressure and the energy exchange densities are comparable to the 
brane energy density Adopting a phenomenological fluid ansatz and generalizations of 
it, we derive a set of exact solutions of the Friedmann equation that exhibit accelerated 
expansion. We find that the effective equation of state parameter for the dark energy 
can exhibit w — — 1 crossing without the presence of exotic matter. 



1 Introduction 



The present cosmological data [I] [2] indicate that a significant percentage of the total 
energy of the universe is attributed to a component referred to as "dark energy". Dark 
energy is the driving force behind the observed accelerated expansion of the universe. 
Its physical origin is largely unknown to this point and the subject of intense theoretical 
speculation and research. Among the numerous cosmological models of dark energy are 
that of a cosmological constant, exotic forms of matter which violate energy conditions, 
such as phantom fields and quintessence as well as modifications of gravitational 
theory itself jS]. A cosmological constant is admittedly the simplest among these attempts, 
although the huge fine-tuning required for its magnitude makes theorists very uneasy [SJ. 
Independently of the challenge posed by the dark energy puzzle, in recent years, the- 
ories of extra spatial dimensions, in which the observed universe is realized as a brane 
embedded in a higher dimensional spacetime (bulk), have received a lot of attention. 
Ordinary matter is confined on the brane but gravity is free to propagate in the entire 
spacetime [7] In these theories the cosmological evolution on the brane is described 
by an effective Friedmann equation [lQj that incorporates non-trivially the effects of the 
bulk. Brane models open up new possibilities for the treatment of standing cosmological 
issues like the observed accelerated expansion [TT]. 

The existence of a higher dimensional embedding space opens up the possibility of 
bulk matter. Bulk matter, since it certainly modifies the cosmological evolution on the 
brane, could possibly be a major contributor to the dark energy. Since, by assumption, 
we cannot probe the bulk, the cosmological evolution on the brane, modified by the 
dark energy through its gravitational attributes, provides us also with a theoretical tool 
to study the bulk. It is also possible that the seemingly phantom-like characteristics of 
dark energy could clXISG clS db result of the way that bulk dark matter affects the brane 
cosmological evolution, even though none of the matter constituents in either the bulk 
or the brane has these exotic features [T^]. 

It turns out that both the bulk pressure and the 05 off-diagonal term of the bulk 
energy-momentum tensor can enter the cosmological equations and have a major effect 
on the cosmological evolution on the brane. Models with energy exchange between the 
brane and the bulk have already been constructed |13j and this kind of interaction seems 
capable of driving the acceleration under certain assumptions. In this article we consider 
the cosmological evolution of the brane in the presence of bulk matter. We adopt a 
phenomenological description of bulk matter in terms of a fluid and derive exact solutions 
of the effective Friedmann equation with an accelerated expansion profile. We find that, 
although the equation of state parameter of the bulk fluid is conventional, for suitable 
values of the bulk energy density and the dark radiation constant, the effective equation 
of state parameter for the dark energy can cross the " — 1" line. 
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2 General Framework 



Consider an Action of the general form 



S = J d 5 xV^G ( 2M 3 R - A + Cf ] ^j + J d 4 xy/=g~ (-a + C 



(m) 
b 



(1) 



Gmn is the 5D metric with signature (— , +, +, +, +) and g^ u the 4D metric on the brane. 
R is the five-dimensional Ricci scalar, while A is the bulk cosmological constant and a the 
positive brane tension. We have included arbitrary matter content both in the bulk and 
on the brane through and respectively. M denotes the five-dimensional Planck 
mass. Varying the action with respect to the metric, we obtain Einstein's equations 

1 1 

Gmn = Rmn — -GmnR = Tmn ■ (2) 

The energy-momentum tensor Tmn resulting from the above Action is of the form 

Tmn = T MN + T M \ — Gmn^ — 9fiu o~ 5{y) ^m^n > (3) 

where T^n results from and T M ^ N results from C^. 
Our metric ansatz will be of the form used in ^1] 

ds 2 = -n 2 (y,t)dt 2 + a 2 (y,t)-f ij dx i dx j + b 2 (y,t)dy 2 (4) 

or 

/ -n 2 (y } t) \ 

Gmn = a 2 (y } t) li3 (5) 

V b 2 (y,t) J 

so that we have a Friedmann-Robertson- Walker brane geometry with a maximally sym- 
metric 3-metric 7^. The fifth dimension is denoted by y and the brane is situated at 
y = 0. Z 2 symmetry is assumed and y takes values in the [—00, +00] line. The met- 
ric functions n(y,t), a(y,t) and b(y,t) are continuous with respect to y but may have 
discontinuous first derivatives at the location of the brane. 

For a bulk energy-momentum tensor T M ^ we shall take the matrix 



-£P 5 P B J \ -n 2 P 5 P B b 2 



-Pb P 5 \ / p B n 2 -n 2 P 5 

= [ P B 8\ , T^ N = P B a 2 1%J ] . (6) 



Note the presence of the off-diagonal component T\ = P5 signifying the flow of energy 
towards (or from) the brane. Note also the anisotropic choice P B 7^ P B , in general. 

For a brane energy-momentum tensor T MN we shall take the matrix 

K»» = S M(~o ) , = 'M ( "0 J ltj ] . (7) 
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Substituting the above ansatze for the metric (0) and for the energy-momentum 
tensor (0), (JZJ) into the equations of motion (J2J), we obtain the equations 




n2 f » <%) / 



4M 3 



a 2 [ a 1 f a! n'\ V ( n' a' \ „ a" n" 



wfTy - - + 2- -- -+2- +2- + 
cr \ a \ a n I b \ n a a n 



a b ( ah 
2_ + _ -2- + 




a b \ an 




AM 3 




A + P B + -f+(p - ^)|7y, (9) 

A, (10) 



a' \ n 2 



4M 3 



{ a' ( a! n'\ b 2 ( a ( a h\ a\ , & 2 1 b 2 , , — \ . . 

where dots denote time derivatives and primes derivatives with respect to y. The curva- 
ture of the maximally symmetric internal space is parametrized by k — — 1, 0, 1. 

At this point we may choose Gauss normal coordinates, such that b(t,y) = 1. In 
addition, we have the freedom to take n(0, t) = no(t) = 1. After theses simplifications 
the junction conditions on the brane, obtained from (|SJ) and Q, are 

a '(+M) = -^(^ + P)> ™'(+M = -^(a-2p-3p), (12) 

where ao(t) = a(0,t) and no(t) = n(0,t). We have used Z 2 symmetry taking a'(— 0, t) = 
-a(+0,t) and n'(-0,t) = -n'(+0,t). 

Considering the last two equations of motion on the brane and using the junction 
conditions, we get 

p + 3(p + p)^ = -2P 5 , (13) 
a 

I + ©' + 4 = (^¥ (ff + rt(2 ^"- 3p) + l^< A ■ (14) 

The first of these equations expresses energy conservation on the brane. For a non- 
vanishing P 5 , energy exchange between the brane and the bulk is present and the energy 
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density p on the brane is obviously not conserved. The second equation will provide 
us with a modified Friedmann equation on the brane. Notice that the "55" component 
of the bulk pressure appears on the right hand side of this equation ((Tj)) affecting the 
cosmological evolution. 



3 Friedmann Equation and Energy Conservation 

We proceed by considering equation (|T4^) which determines the time evolution of the 
scale factor on the brane. We shall assume an equation of state p = wp to hold between 
the energy density and pressure of matter on the brane. Dropping the subscript "o" and 
defining (3 = (24M 3 )" 2 and 7 = a/3, we get [T3] 

a . / d\ 2 k „ s „ 2 Pb A 



where 



<T 2 



We shall make the assumption that the effective 4D cosmological constant A vanishes 
when no matter is present. Thus, we shall proceed taking A = 0. We shall also be 
interested in the low density case, i.e. 

(l + 3w)p 2 « (1 -3w)pa. (17) 

Our second order evolution equation ()15|) . with A = 0, can be cast in the form of two 
first order equations, one of which is analogous to the Friedmann equation of standard 
cosmology. This can be achieved by the introduction of a dark energy variable x(t). Our 
new set of equations is 

a \ 2 k Pr /x 

+ - = 2 1P + (3p 2 + x (18) 



a l a 



12M 3 



a ( Pb \ Pb 



^ + 4 a^-24^J-lW- 4P ^ + ^^°- (19) 

The pair of these equations 1 is equivalent (JT5J). 

In addition to these two equations we also have the energy conservation equation 
(fT3|). Assuming the low density approximation, our set of equations is 

p + 3p- (l + w) = -2P 5 , (20) 
a 



1 Comparing the first of these equations to our original evolution equation (|f 5|) . we have 



-X- (3w + 2)f3p z -(3u> + l) 

(X 



Differentiating l|f 8() and using this equation proves (|f 9|) . 
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f a\ 2 k Pr , x 

UJ + ^ = 2w + ^-i2i§5' (21) 

* + 4 S = I^ +4P57 ' <22> 

Equation is the Friedmann equation describing cosmological evolution on the 
brane. The auxiliary field x incorporates non-standard contributions (dark energy) dif- 
ferent than the regular matter content of the brane. The bulk matter contributes to the 
energy content of the brane through the bulk pressure terms (P B = T 5 ( f } ) that appear 
in the right hand side of the Friedmann equation. In addition to that, the bulk matter 
contributes to the energy conservation equation (J20|) through P 5 = T° that accounts for 
energy moving from the brane to the bulk or the opposite. The functions Pb, P5 are 
functions of time corresponding to the values of Pb{u, t) and P$(y,t) on the brane. The 
energy-momentum conservation VmT% = cannot fully determine Pb and P5 and a 
particular model of the bulk matter is required. In what follows we shall model the bulk 
matter in terms of a phenomenological regular bulk fluid ansatz that leads to analytic 
solutions. 



4 Digression on Bulk Matter 

In this section we shall derive the conservation equations for the components of the bulk 
energy-momentum tensor and investigate the constraints of the form of P5 and ?gon the 
brane for various ansatze and in the particular case of a fluid. Consider the conservation 
of the energy-momentum tensor 

VrpM ri , Q rpM I pM rpR -pj? rpM ri ^OQ"! 

M-L N - U °M 1 N + 1 MR 1 N ~ 1 MN R ~ U • V Z6 ) 

The "0" component of O is 



Pb + 3-(p fl + Pb) + 3n 2 P 5 (- + -)+ n 2 P^ + 5(y) ip + 3-(p + p) ) = . (24) 
a \n a J [a J 

Integrating around y — and using the Z 2 symmetry (P 5 (+0, t) = — P 5 (— 0,t) = Ps(t)), 
we reobtain the "05" equation of motion on the brane 

p + 3-{p + p) = -2P 5 , 
a 

and, in the bulk, the equation 

Pb + 3-(p B + P B ) + 3n 2 P 5 (- + -)+ n 2 P' 5 = . (25) 
a \ n a J 



Similarly, the "5" component of (j23j) gives 

P' B + P 5 + P 5 ( - + 3- ) + - (Pb + Pb) = . (26) 



n a n 



In the limit y — > +0, these equations become 

h + 3^(p B + P B ) + P' b = ^{2a - p - 3p)P 5 , (27) 

P 5 + 3^P 5 + T' B = ^^(a -2p- 3p)(p B + P B ). (28) 

These two equations as they stand cannot limit the freedom of choice for P B (t) and 
Ps(t), since the unknown functions P 5 ' and P B appear in them. Nevertheless, specific 
models for the bulk matter will greatly eliminate this freedom. 

As a concrete model let us consider the energy-momentum tensor for a classical fluid 
in motion, moving at a speed v along the fifth dimension 

T% = P B 0) 5% + (pf + P<°>) U M U N , (29) 

with 

U° = (1 - v 2 )~ 1/2 , IP = 0, U 5 = v(l - v 2 )- 1 ' 2 . 
Assuming a fluid-type equation of state 

P® = upf , (30) 

we get 

/ u-n 2 (l + uj) v(l+u) \ 
T M N « pg } u8) + 0(t; 2 ) . (31) 

\ -n 2 w(l + u;) u + v(l + u) J 

In terms of the standard variables of the previous sections, we have 

T° = —p B = -pf [n 2 (l , T° 5 = P 5 = pg> v (l + J) , (32) 

T 5 = -n 2 P 5 = -n\{\ + ui)pf , T) = 6)P B = upftfj , (33) 
T\ = P B = pf[u + v(l + u)} . (34) 



Substituting these quantities into equations (}2*7j) and (J2EJ) we obtain two equations in- 
volving p B and . Cancelling the spatial derivative among them, we obtain a single 
first order equation for p B (0,t) 7 namely 

pf [u + v(l + u)} + pf |3^(1 + u) + 3~u(l + w) 2 

+ ^[3^(l + ^)(p + 3p-2a) + (l+u) 2 (a-2p-3p)]} = 0. (35) 
The lowest approximation solution to this is 

Pb « ^b) ( 1 + 0(w/M 3 ) + 0(v 2 ) ) . (36) 



This corresponds to 



Pb = C 1 + •-), PB = u-Jj^ (1 + -..), (37) 



P 5 = W (l + W )^(l + 0(, 2 )) , P B = + -..) • (38) 

The dots stand for 0(va/M 3 ) and 0(v 2 ). 

From a phenomenological point of view it is possible to generalize this ansatz depart- 
ing from its fluid interpretation. For instance, inserting time-dependence in the parameter 
v would introduce an extra v term in (|33j) but the lowest approximation (j3*Uj) would not 
be affected. 



5 Exact Solutions and Late Time Acceleration 

Let us now return to our set of cosmological evolution equations in the low-density 
approximation (3p « 7 

^ 4 K x -^)^ 4/3 ( p+ i) n+ i^ K47n+ i^' (4o) 



p + 3-p(l + w) = -2P 5 . (41) 
a 



Using the associated equation 



a 



a 



(3w + 2)Pp 2 -(3w + l) 1P - X « -(3w + l) 7 p-x, (42) 



we can write a general expression for the deceleration parameter q, which for k = is 



a\ 2 fa\ x + (1 + 3w) r yp 



q = - - - = - v L ir . (43) 

W W x + 2 7 p - 

Accelerating behaviour emerges when g < 0. 

The time dependence of the scale factor a(t) is determined by the Friedmann equation 
(|39j) which, ultimately, will have a general form 

/ • \ 2 at n 

u) = c ° + E c » an + E c -« a_n • ( 44 ) 

^ ' n=l n=l 

Assuming that we have an expanding behaviour, at late times the right hand side will 
be dominated by one term, namely 

2 

CX . (45) 



Owl 



The resulting late time scale factor is either 



a(t) = a(to) [l - V -\fC v W*)? lv (t-t Q ) 



for v < 0, or 



a(t) = a(t ) 



t r — to t 

for v > 0. The time t r is the "big rip" time for which the scale factor blows up. Accelerated 
expansion can be achieved in a trivial fashion when the dominant term is a constant. 
In general, accelerated expansion requires q = — 1 — I//2 < . It is interesting to see 
whether bulk matter can induce accelerated expansion in the case that no cosmological 
term is present (A = 0) and ordinary matter is in a standard w = or w = 1/3 phase. 

The set of cosmological equations (|39|). ()4L)j). (jTTj) can be solved exactly for the 
following general ansatz for the bulk matter pressure on the brane (see also [T2] ) 

P B = Da" , P 5 = ^ (46) 



We shall proceed to obtain the resulting exact solution but we will justify this ansatz 
later. Substituting (|4l)|l in equation (JJUJ), we obtain 



X + A- a (x- 5a» - P -5a v - 7 = 



or 

a 4 (z/ + 4) (fj, + 4) 

where we have defined 5 = D/2AM 3 . Similarly, we can proceed with the integration of 
(jUJ and get 

C 2F 

p = _ - _ - a » . (48) 

P a 3 ^"-) [3(1 + w) + n] K J 

Substituting (|4*7jl and (}4*SJ) into the Friedmann equation (|39p. we can write it in a con- 
ventional form as 

+ ~2 = -Z-G* Peff , (49) 
CI o 

where = 37/47T = 3cr/47r(24M 3 ) 2 is the 4.D Newton's constant and the effective 
energy density p e ff stands for 

C C/2 7 _ 25 v 2(3w- l)F ^ 

Peff ~ a 3(l +W ) + fl 4 l{y + 4) fl + (/i + 4) [ 3(1 + w) + ^] • ^ 

Notice that for w = 1/3 the T\ decouples from p e //- 
The acceleration behaviour can be deduced from 

o 7 (l + 3w)C g „, (y + 2) „ 2F 7 (l-3 W )(^ + 2) „ 
a a 3 ^-) a 4 (i/ + 4) (// + 4) [3(1 + w) + /x] 



S 



For v and \x greater than the standard matter exponent —3(1 + w), the bulk generated 
terms will dominate at late times. The corresponding late time deceleration parameter 
comes out to be 

' v>n -1 - f 



Q 



fi = v -1 - ~ • (52) 



Note that the late time deceleration parameter comes out independent of the particular 
ratio of the ansatz parameters D/F. This means that even if one of the two bulk pressure 
terms vanishes, the late time deceleration behaviour would be the same. This is not true 
of course for p e ff which is proportional to —45 + ■ 

A justification for the ansatz (J46|) that we employed can be found in the framework 
of the simple fluid model that we analyzed in the previous section. There, we found that 
an approximate solution for a conserved bulk energy-momentum tensor lead us to (|37p 
and PS]1. The second of them reads 

P s = v(l+u)j^, P B = ' 

As we mentioned there, we can adopt a phenomenological point of view and consider a 
velocity parameter that depends on time. A measure of temporal development on the 
brane is given by the Hubble parameter H — a/ a, so it would be natural to take v = (H, 
where ( is a phenomenological parameter. With this choice, the above approximate 
solution takes the form 

*»=^. « = C(l+<*)^(i)- (») 
Thus, for our phenomenological bulk fluid we obtain the effective energy density 2 

C C/2 1 2cu C B 2C(3w-1)(1 + ^) C b 

Peff a 3(i+ w ) + a 4 7 (1 - 3cu) a 3 (!+-) (1 - 3w) [3(w - w)] a 3 (!+-) 1 ' 

and 

a 7(l + 3w)C C (H-3w) C B 2Ct(1 + w) (1 - 3w) (1 + 3w) C B 

7 + 2^-- — : ^7TT~^\ + 



a a 3 ( 1+ ^) a 4 (1 - 3u) a 3 ^) (1 - 3u) [3{w - u)} a 3 ^) ' 

(55) 

For u) < 0, the bulk term is dominant at late times and gives a decceleration parameter 

» = 5 + 3|. (56) 
The corresponding acceleration range is 

q<0^-l<uJ<-\ (57) 



2 Where C B = and C = 24M 3 C_ 
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Both contributions to the Friedmann equation different than the brane energy- 
density, namely, the dark radiation term C/a 4 and the bulk terms, can be interpreted as 
dark energy. The effective equation of state parameter for the dark energy w^J can be 
computed from the prescription [ 



w 



(D) 
eff 



^ ld\n(5H 2 ) 
3 din a 



(58) 



where 5H 2 = H 2 /Hq — Q m a~ 3 accounts for all terms in the Friedmann equation not 
related to the brane matter, for which we have taken w = 0. In our case, equation 
becomes 



w 



(D) 
eff 



-1 + 



1^+3^ -4c+fC 7 



C B (l + z 



,3w-l 



c + 



1 



(l-3w) 



-4c + f(7 



C B (l + z 



)3w-l | 



(59) 



The denominator is proportional to the effective dark energy density and should be 
positive. In order to achieve w^J < —1, we need 



4C +3- 



c 



4- (1+c) 

-4c + -(^ >- 

3 c 



C B (l + z) 



3a;-l 



< 0. 



(60) 



;i - 3c) 

Taking into account that the denominator is positive, we obtain the following lower 
bound for the above expression: 



AC +3- 



4- (1+w) 

-4c + -C7^ 

3 uj 



C B {1 + zf"- 1 > 



1 -3uj)C. 



(61) 



;i - 3c) l 

We want the lower bound to be negative. As already mentioned, in order to get acceler- 
ated expansion, oj must be within the range — 1 < uj < — |. For this range of values, we 
obtain 1 — 3c > 0. Thus, in order to get a negative lower bound, we must assume C < 0. 
Clearly, the existence of such a lower bound is a necessary but not sufficient condition 
for w = —1 crossing to occur. However, with proper values for the parameters ( and C B 
of the bulk fluid the left hand side of ()61|) can assume negative values as the redshift z 
varies. 

Summarizing, we presented a treatment of the cosmological evolution of braneworlds 
for general bulk matter content. In our framework the bulk matter pressure and the 
energy exchange between bulk and brane are comparable to the brane energy density. 
We adopted a phenomenological description of bulk matter in terms of a fluid and using 
the ansatz T B 5 ~ a we derived exact solutions of the Friedmann 

equation that exhibit accelerated expansion. We found that the effective equation of state 
parameter for the dark energy can exhibit w = —1 crossing without the use of matter 
that violates the energy conditions. 
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